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Abstract We report on progress concerning the partial wave analysis of higher cor- 
relation functions in conformal quantum field theory. Q 



1 Introduction 

Partial wave analysis (PWA) is a powerful tool in conformal quantum field theory. 
It gives not only information about the field content and the operator product expan- 
sion (OPE) of a model EG), but can also be used for probing the positivity of the 
inner product induced by the correlation functions (Wightman positivity) ifTOl . 

Positivity is difficult to establish because it is a nonlinear property. It also neces- 
sarily involves correlation functions of any number of fields iFTol . The most promi- 
nent example is the classification of central charges below 1 of the Virasoro algebra. 
An example in four spacetime dimensions (4D) is the result that conformal scalar 
fields with global conformal invariance (GCI, |[T3ll ) are necessarily Wick squares of 
free fields 111 111 , and cannot couple in a nontrivial manner to other fields J2). 

While conformal PWA for 4-point functions is well understood Q, we intend to 
develop methods for higher correlation functions. The basic task is to decompose a 
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correlation function of conformally covariant fields into a sum over partial waves 

{n,(j>i(xi)...<j) n (x n )n) = Y,(fi,fr(xi)---Qk-i(xk-i)nxfa(xk)---Qn(xn)fi), (i.i) 

X 

where Tlx m me projection to the subspace of the Hilbert space which carries the 
irreducible representation X of the conformal group. A projection can be inserted in 
any position within the correlation, so that the n-point partial waves depend on « — 1 
representations, where the first and last projections are redundant because they are 
fixed by the first and the last field. 

In principle, the non-vanishing partial waves give information about the contribu- 
tions to the OPE of two or more fields (8j. Since a projection is a positive operator, 
each partial wave contribution of the form 

(xi) .. . 0(x„)n A 0(x„ +1 ) . . . f (a&,)fl) 

must separately satisfy Wightman positivity (i.e., after smearing with test functions 
f(x n , . . . ,xi)f(x n +\,. . . ,X2n) it must yield a non-negative number which is the norm 
square of the vector 17(0 <E> ■ ■ ■ ® (j)') (f)Q). More generally, partial waves are subject 
to Cauchy-Schwartz type inequalities. 

Now, partial waves are to a large extent determined by conformal symmetry, be- 
ing solutions to eigenvalue equations for the Casimir operators of the conformal 
group. Therefore, the positivity requirement reduces to the positivity of a numeri- 
cal coefficient, the partial wave amplitude, which multiplies a model-independent 
partial wave function ifTOl . 

Conformal PWA is by now mostly limited to 4-point functions, because the 
higher partial waves are not sufficiently well known. Even for 4 points, the de- 
termination of partial waves in 4D required a considerable effort (6j. Moreover, 
the decomposition of a given correlation function into a known system of partial 
waves may not be a straight-forward task without a suitable notion of orthogonality 
between the partial waves. Some progress was made in iflOl giving a systematic ex- 
pansion formula for scalar 4-point partial waves, and in |[T4l for a suitable notion of 
orthogonality. 

In this note, we report some further intermediate progress. In Sect.|2] we present a 
power series representation Eq. ( 12.4b for general n-point partial waves in two space- 
time dimensions (2D) for all n, extending known formulae for n < 4. In 4D, how- 
ever, such an expansion seems unrealistic because of the complicated structure of 
the higher-order Casimir operators which the partial waves must diagonalize, and 
because the partial waves are no longer unique. 

In Sect. [3] we therefore present an alternative to the actual decomposition Eq. 
( 11.11 ). which is applicable also in 4D. The idea is a successive reduction of n-point 
functions to n — 1 -functions, in terms of local linear maps (j>\ (xi)^(x2)i2 — > §x (x)Q 
selecting each contribution to the OPE of the last two (or the first two) fields in the 
correlation. Our main result is the characterization of these linear maps as partial 
differential operators that intertwine the respective representations of the conformal 
group. This property is encoded in Eq. ( 13.4b . which is subsequently solved. Acting 
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on the correlation functions, the intertwiners effectuate the desired reduction. As we 
shall see, this method is applicable only for representations of integer scaling di- 
mension (otherwise, the differential operators would have to be replaced by integral 
kernels Q, and locality would become a nontrivial issue). 

This method is therefore well-suited for QFT with global conformal invariance, 
where all correlation functions are rational functions |[T3l , We shall apply it in Sect. 
|4]to address the problem of positivity of a class of "exotic" higher (n > 6) correlation 
structures of twist 2. The motivation is the following. 

Twist-2 contributions in free field theories above the unitarity bound arise from 
quadratic Wick products such as : (p*(x{)(p{x2) : orxi2^ l '.y/{xi)Y >i y/{x2)' , in which 
each factor can be contracted "only once", so that both variables can only have 
poles w.r.t. one other variable. In contrast, the exotic structures contain so-called 
double poles, thus indicating a nontrivial theory. These are strongly constrained by 
the conservation laws for twist-2 fields ifTTl . allowing for a classification ||4]. In 
particular, they cannot arise in correlations of less than six fields. While the exotic 
structures satisfy all linear properties, it remains an open problem whether they are 
compatible with positivity. 

First steps of the positivity analysis of the simplest exotic structure will be re- 
ported in Sect. |4] 



2 Higher chiral partial waves 

Irreducible representations X of the conformal group are eigenspaces of the Casimir 
operators. Thus, correlation functions with projections onto irreducible subrepre- 
sentations inserted: 

(a , </>! fa )n Al fc(* 2 ) • • • n Vl ^{xi)n Xi ■ ■ ■ ^(x^n^u^a) (2.1) 

are eigenvectors of the corresponding differential operators arising by commuting 
the conformal generators with the fields. Partial waves are, by definition, solutions 
to the same eigenvalue differential equations, with some standard normalization. 
These are "universal" in the sense that they are completely determined by conformal 
symmetry. They depend on the sequence of representations (jU/)j=i,.. n of the fields 0,- 
in Eq. ( 12.11 ). and on the sequence of representations (A/);=i,..„_i of the projections, 
where Ai = ~p l and A„_ i = ji n are redundant. 

The projected correlations Eq. ( 12.1b are multiples of the partial waves. The coef- 
ficients contain model-specific information, and Wightman positivity can be formu- 
lated as a system of numerical inequalities on the partial wave coefficients 1 10']. 

The conformal Lie algebra in 4D, so(4, 2), has three Casimir operators (quadratic, 
cubic and quartic in the generators). In contrast, the conformal Lie algebra in 2D fac- 
torizes: so(2,2) ~ s/(2,R) ®sl(2, K), and each s/(2,R) has one quadratic Casimir 
operator. For this reason, the Casimir eigenvalue differential equations are much 
simpler (both, to write down and to solve) in 2D. 
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The relevant positive-energy representations of sl(2,M.) ©s/(2,R) are parame- 
terized by the chiral scaling dimensions d±, such that d + + d_ is the total scaling 
dimension, and d + — d- the helicity. 

Because of the chiral factorization of the conformal group, also the partial waves 
factorize. In the sequel, we display only chiral partial waves as functions of either 
x+ = t+x orx- = t — x, and suppress the subscript. Thus, a general projected cor- 
relation function in 2D has the form of a product of two chiral functions 

(a , 0i (xi )n ai $ 2 ( X2 ) ■ ■ ■ n ai _, ■ ■ ■ 0„_ , (*„_ , )n Un _ l 0„ ( x „)n) (2.2) 

where the chiral fields 0,- have chiral dimensions dj, and TI a are the projections onto 
the chiral representations with chiral scaling dimension a. In particular, a\ = d\ and 
a n -\ — d„ are fixed. 

The Casimir eigenvalue equation for the projector insertion Tl aj reads 

£ *jk d J d k + 2 E d /(*M) + £ d * _ (£ d *) 

i<j<k i<j,k i<k i<k 

= { a i - a l){- ■ ■ naA+l • ••}, 

which is equivalent by conformal invariance to 

( £ xjJA + 2 E d i +£4-(£4)V-<fc (xdn ai ...) 

j<k<i j,k<i k<i k<i 

= (ai-af)(...<l>i(xi)n a . ...). 

In principle, these equations can be reformulated in terms of n — 3 independent 
conformal cross ratios. It turns out convenient to choose 

Xkk+l x k+2k+3 

Uk = ■ 

Xkk+2*k+lk+3 

We have worked out the invariant differential equations for n < 6 points: Let 
*2j2l = f(ui,u 2 ,u 3 ) 

'— ' di+d 2 -d3 di+d 3 -d 2 d 2 +d 3 -di A+d 5 -d b d 4 +d 6 -d 5 d s +d 6 -d A ' 
A 12 x 13 A 2 3 ' x 45 A 46 x 56 

Then (with the Euler operators Ej = uid Uj ) 

(E l +d 3 -a 2 )(E l +d 3 +a 2 - \)f = u\(E] + E 2 )(E { + d { -d 2 + di)f, 

(E 2 -a 3 )(E 2 +a 3 -l)f = u 2 (E 2 +E 1 )(E 2 +E 3 )f, 
(Ei +a?4 -04) (£3 +<i4 + «4 - 1)/ = Ui(Ei +E 2 ) (Ei +de-d$ +d^)f. (2.3) 

(The cases n < 6 are covered by admitting the trivial field 1 of dimension 0.) This 
system is obviously symmetric under hermitean conjugation 1,2, ... ,6 —¥ 6,5,..., 1, 
It can be recursively solved as a power series with leading powers u" 2 di u£u^ dA . 
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From the solution with n < 6, we have extrapolated the general power series 
expansion for all n, as follows. By default, we put oq = a n := 0, and (.q = £„_2 '■= 0. 
Proposition 1: The general chiral n-point partial wave is 



y m^lThl^lZ^h^i ff ^ (24) 

This formula has a remarkable "short-range" feature: It involves only coordinate 
distances xij with j = i + 1 or i + 2. The powers of andx,, + 2 depend only on the 
dimensions of the fields 0,-, respectively 0;+i, and their adjacent projections, 
apart from the summation indices I. The same is true for the numerical coefficients. 

For n = 3 points, this is just the 3-point function. For n = 4,5,6 points, we have 
derived this formula by solving the differential equations Eq. ( 12.31 ) for the Casimir 
eigenvalues. For n = 4, the sum is a hypergeometric series, and Eq. ( 12.41 i coincides 
with well-known formulas. 

One way to prove the Eq. ( 12.4b for all « is an application of the method discussed 
in the next section. There, we introduce "intertwining" differential operators 1 o 
E/ t with the distinguishing property that they annihilate all partial waves carrying 
the "wrong" representation a ^ h, and reduce the n-point partial wave carrying the 
representation a = h to an (« — l)-point partial wave with the first pair of fields 
replaced by 0o of dimension h. 

Therefore, it is sufficient to show that this is true for our "candidate" partial waves 
Eq. ( 12.41 ). With Eq. ( 13.181 ), we have to apply the differential operator 

d 1 +d 2 _f V- (q-b) p (p + b) q - p - . \ di+d 2 



E h = E h o Xl \ 2 = ( L — ~, ~, ^(-d 2 )' 

p+q=h P- C l- 



ox 
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where b = d\ — di, to Eq. ( 12.4b . and then equate x\ = X2- The result must be 
times the reduced partial wave. 

To do this, we have to exhibit all terms that involve x\ or x 2 . Eq. (12.4b can be 
arranged as xj~^ 1_< * 2 times the sum Yd 2 ,—,t -3 over 

xn Y f x 2 3\ b f X23Y (a + b) e (a + c)eu e . . 
— — L ink ^ — x remaining factors, (2.5) 

X13X24J \xi3 J \x24J t> V\2a)t 

where a = a\ , b = d\ — dj,, c = 03 — dj, + £2- Notice that for each 1%, the sum over I 
is a 4-point partial wave where the 4 th field has dimension 03 + £3. Thus, knowing 
that Eq. (12.4b correctly reproduces the 4-point partial waves, and that 1 o Ef, reduces 
4-point partial waves to 3-point partial waves, the same must be true for the higher 
partial waves. 

However, we have not been able to evaluate the result of E^ 1 ' 1 ' 2 on the power 
series Eq. ( 12.5b . and verify this indirect conclusion by a direct computation. Only 
for n = 3 this can be done by the following argument. For « = 3, one has c = in 
Eq. ( 12.5b , only £ = contributes, and there are no "remaining factors". Then 
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(i) Because E' h [ '' 2 is a differential operator of order h, it annihilates the 3-point 
function whenever h < a, due to the surviving factors of x\2- 

(ii) Writing Xu = — — and performing the binomial expansion of its pow- 

ers, E^' dl can easily be applied. It is then seen by inspection that the resulting series 
is symmetric under the exchange a <-> h. Therefore, it also vanishes whenever h > a. 

(iii) When h = a, all derivatives must hit the factor x" n - That the result is the 
2-point function, is then obvious. 

For n > 3, the Leibniz rule produces multiples sums which are not easy to handle. 
But a trick helps: The sum in Eq. ( 12.5b equals 2F1 (« + b, a + c;2a;u). We then use 
the identity 

2a— 1 / \ b / > 

34 I x 23 \ I *23 1 „ , . , „ \ 

1 — I — I ■2Fi(a + b,a + c;2a\u) 



r(2a) 



r(a + c)r(a-c) 



P dx{x x -x)-"- h {x 2 -x)- a+b (x 3 -x) a+c - l {x-x 4 )"- c -\ 

J XT. 



which can be established by direct computation: namely, the change of variables t = 
*24fe-*) yj e ids precisely the standard integral representation of the hypergeometric 

.1:34 [X2 X) 

function. 

Therefore, each term Eq. (12. 5t is, as far as its dependence on x\ and x 2 is con- 
cerned, an integral over a 3-point function. Thus, we only have to evaluate E d> ' d2 on 
a 3-point function, which can be done as before. The remaining integral is again of 
the hypergeometric type (after the change of variables t = — ^ ). and reproduces 
precisely the necessary "leading" factors for the (n — 1) -point partial wave Eq. ( 12.4b - 

From this, we conclude that Eq. ( 12. 4t indeed is the correct power series expansion 
of general n-point chiral partial waves. 



3 Intertwining differential operators 

Let 0i and 02 be two conformal fields transforming in representations [i\ and /io. 
We shall determine differential operators w.r.t. x\ and x 2 such that 

fa(x) :=l x °Exfa(xi)fo(x2) (3.1) 

transforms like a conformal field in the representation X. Here, l x is the evaluation 
map i x (f) = f(xi ,x 2 )\ xi =x 2 =x = f(x,x). 

It will become clear below that such operators exist only when the scaling di- 
mensions satisfy dx—d\—d%^L 7L. They can therefore be expected to be exhaustive 
(w.r.t. A) only in a globally conformal invariant (GCI) theory. 

Such operators have been presented previously J9] Sect. VLB] for the special case 
of 0i and 02 being two (complex conjugate) canonical scalar massless Klein-Gordon 
fields of dimension 1, in order to extract the current, the stress-energy tensor and 
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higher conserved symmetric traceless tensor fields from : 0*0 : . The same operators 
actually can be used also for scalar biharmonic bifields V{x\,X2) which collect the 
twist-2 contribution in any product of two scalar fields of equal dimension ifTTI . 
where biharmonicity, i.e., the wave equation w.r.t. both arguments is exploited in an 
essential way. We shall reproduce these operators, but there will be additional terms 
including the wave operators, so that Eq. (13. U is true without using the equation of 
motion, or biharmonicity. 

By conformal covariance, the assumed transformation behaviour of Eq. (13. U im- 
plies 

h ° £a (0m b>)& , 01 (*i ) 02 (xi)Q ) = <5 AAt (<pn {y)Q , A (x)Q), 

i.e., the operator annihilates all 3-point functions with fields in the "wrong" repre- 
sentation. In particular, if applied to the vacuum operator product expansion (8) 

fr(jci)fe(x 2 )fl L I "dxKf! m (x 1 ,x 2 ;x)(l> ll (x)Q, 
n J 

where K^ >fi2 are certain integral kernels, it will annihilate all contributions fi^X, 
and if applied to a correlation function, it will annihilate all partial waves with /I ^ A 
in the 1-2-channel, and reduce the contribution with yt = A to an n — 1 -point partial 
wave. Thanks to the latter feature, one can perform a partial wave analysis without 
actually knowing the partial waves, cf. Sect. [4] 

Let us now proceed to determine the differential operators. 

For definiteness, we specialize to ji\ = jj.2 to be scalar representations of dimen- 
sion d\ = d.2 = d. In this case, only symmetric traceless tensor representations A 
can occur (8j. It is convenient to write A = (k,L) where L is the tensor rank, and 
2k the "twist", such that the scaling dimension is d = 2k + L. The unitarity bound 
requires K > for L = 0, and K > 1 for L > 0. We write a symmetric traceless 
tensor as T(v) = T^-^v^ ...Vn, which is a homogeneous polynomial of degree 
L in the polarization vector v. Tracelessness is equivalent to the harmonic equation 
□ v r(v) = 0. Eq. (13.1b implies that E K i is a harmonic homogeneous polynomial of 
degree L in the polarization vector v. The harmonic part of any polynomial in v is 
uniquely determined O, so it is sufficient to know E K i up to terms involving v 2 . 

Let T = P tll D, M|, v ,^i be the generators of translations, dilations, Lorentz and 
special conformal transformations, respectively, and 

the commutation relations with covariant ("quasiprimary") fields, where t kL = d 
for the translations, = (xd + dx ) for the scale transformations, = xAd + vAd v for 
the Lorentz transformations, and = 2x(xd) — x 2 d +2(v(xd v ) — (xv)d v ) +2d^x for 
the special conformal transformations. For the tensor representations, dx =2k + L, 
while for the scalar representations = ji2 the v-terms are absent and d^ = d. 
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Commuting the generators with Eq. ( 13.11 ), the assumption that <j> x transforms in 
the representation A is equivalent to the intertwining relations 

ioE x o{t^+t^)=t^oioE x . (3.2) 

In the case at hand, we make an ansatz 

E x = E KL (x i ,d i ,v)o(x 2 n ) d . (3.3) 

Notice that by virtue of the pole bounds |13], any correlation function of 0i (xi)(j) 2 (x 2 ) 
is not more singular than (x\ 2 )~ a \ so that the differential operators E k l act on a reg- 
ular function, and the subsequent evaluation i x is possible (provided E K i is regular). 

Next, we evaluate the intertwining relations Eq. ( 13. 2\ . They tell us in turn: 

Translations: (d\ + d 2 )E K i = 0. Thus the differential operators do not involve 
the coordinate x\ +x 2 . Since E K i is followed by the evaluation map l v , we may also 
assume that it does not involve the difference coordinate x\ — x 2 , hence E k l involves 
only derivatives and the polarization vector v. Let us denote by V, the derivatives 
with respect to the "variables" <9 ( - of E K i(d\ ,d2,v). 

Scale transformations: (<?iVi + d 2 V 2 )E K L = (2k + L)E k l. Thus, E K i is homoge- 
neous of degree 2k + L in the derivatives <?,-. 

Lorentz transformations: (<9i A V \ + ^ A V2 + v A d v )E K i = 0. Thus, E k l is a 
Lorentz scalar. It is therefore a function of (<3,<3/), (v<9,) and v 2 . Together with the 
known homogeneities in v and in <?,-, it can be a polynomial in the derivatives only if 
K is an integer. This is in perfect agreement with GCI because tensor-scalar-scalar 
3-point functions are rational only if the twist 2k is even. 

Special conformal transformations: While the previous intertwining conditions 
gave information about the gross structure of E K i, the special conformal transfor- 
mations yield a differential equation that specifies the operators completely. 

Proposition 2: Given the previous specifications ofE K i(di ,d2,v) in Eq. ( Ii.il ) as 
homogeneous polynomials (of degrees depending on the parameters K and L), the 
intertwining condition Eq. ( Ii.2| ) is equivalent to 

(2(5 1 V 1 )V 1 - d x V\ +2(d 2 V 2 )V 2 - d 2 vt)E KL {d u d 2 ,v) = 0. (3.4) 

One may directly solve these equations with a polynomial ansatz for Ekl with 
the specified homogeneities. A more systematic way is to write 

E KL (di ,d 2 ,v) = fab)* ■ [((yd! ) + (vd 2 )) L ■ e KL (p,q, r)] Q (3.5) 

where p = q = ^-y, and r = [^|~^|] . Clearly, e Kh must be a polynomial 

of degree at most L in r, and degree of at most K in p and q. The notation [P(v)]o 
stands for the harmonic part of the polynomial P(v). The variable v 2 does not appear 
explicitly, because the harmonic part [v 2 <2(v)]o = for any polynomial Q J3|. 
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With this ansatz, the differential equation Eq. ( 13. 4t turns into the system of three 
PDE for e KL (p,q,r): 

(L(L - 1) + (1 - r 2 )d 2 + 2k(L - rd r ) + 2{pd p - qd q )d r ) e KL = 0, (3.6) 

[4(pd p -l)d p - q(K- pd p - qd q )(K- 1 - pd p - qd q ) + 

+2(K- pd p - qd q )(K- 1 - pd p + qd q + (r- l)d r )] e KL = 0, (3.7) 

[4 (qd q -l)d q -p(K- pd p -qd q ){K-\- pd p - qd q ) + 

+2{K-pd p -qdq)(K-l+pd p -qdq + (r+l)d r )] e KL = 0. (3.8) 

One may repeat the same strategy in 2D. In this case, the intertwining operators 
factorize into two chiral operators, labelled by the chiral dimensions h±. These are 
polynomial functions in the chiral (one-dimensional) partial derivatives <9i and d2 ■ 
Following the same line of arguments as in 4D, one finds the chiral intertwining 
condition 

(d 1 V 2 l + d 2 V 2 2 )E h (d 1 ,d 2 ) = 0, 

where E/,(di,d2) is a homogeneous polynomial of degree h. Writing E/, = (di + 
di) h ■ e/i( gjrpg: )» this reduces to the differential equation for e/,(r) 

(h(h-\) + (\-r 2 )d 2 )e h {r)=Q 1 

which is exactly the same as the case K = 0, L = h of Eq. ( 13.6b . 

Notice that in 4D, representations (0,L) with L ^ are below the unitarity bound. 
Such representations must not contribute to a correlation function. Thus, any admis- 
sible correlation function must be annihilated by the operators 1 o Eql. The solution 
for K = is 

e 0L (r) = (l-r 2 )d r P L - 1 (r), 
where P n are the Legendre polynomials. Using Eq. ( 13.51 ). this gives 

(3.9) 



E 0L (d l ,d 2 ,v)= £ ^k^k\{yd x )P(- V d 2 f 



or (in the chiral case) 



{l)p {P)q ^p. 



E h (dud 2 )= I ™f^Ld?(-di) q - (3.10) 

For K > 0, we may expand 

e KL (p,q,r)= £ P" V ' e KL-mn{r) ■ 

m,n>0,m+n<K 

Then Eq. (I3.61 l must hold for each term p'"q n e K L z mn{r) separately, giving 

(( 1 - r 2 )d 2 - 2Krd r + 2{m - n)d r + L(L + 2k- \))e KL , mn (r) = 0. (3.11) 
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This equation involves only the difference m — n=: 8. It is solved by polynomials 
of degree L with the symmetry f KL . s (r) = (- 1 ) L f KL ;-s{-r)\ 

ftLfi W = (K - S) L ■ 2 Fi ( -L,L + 2K- 1; K- 8; ^) . (3.12) 
Thus, to solve Eq. ( 13. lit it remains to determine only the coefficients in 

e KL;mn { r ) = C KL\mn ' fKL;m—n { r ) ■ (3.13) 

Indeed, the remaining Eq. ( 13.7b and Eq. ( 13.8b turn into the recursive system 

4(m 2 - 1 )c KL - m+Ln +2(K-m-n)(L+K-l-m + n)c KL - mA 

-(K-m-n){K-m-n+l)c KL -j n ,n-\ = Q, (3.14) 

4(« 2 - l)c KL ;m,n+i+2(K-m-n)(L+K- 1 +m-n)c KL -m,n 

—(K — m — n)(K—m — n+ l)c KL - m ^\^ n = 0. (3.15) 

Here, we have used the fact (|1 Eqs. 15.2.14 and 15.2.16]) that the differential 
operators 

± _ {r^\)d r +K-\T8 
KL > 5 ' L+k-It8 
act as raising and lowering operators for the parameter 8: 

A kL,S f<cL;8 = /k£;5±1- (3.16) 

We conclude: 

Proposition 3: The intertwining differential operators in Eq. rti.il ) are given by 



E KL = £ c KLtrnn {d 1 & 2 ) K - m - n D^n" 2 {vd l+ vd 2 ) L f KL .^ n ( V ^ ; "I 2 ) °( x n) d 

m+n<K L Vdx+Vd2 Jo 

where [. . .]o stands for the harmonic part with respect to v £ R 1,3 , f/ze polynomials 
ftcL;m-n ore given by Eq. A3.12\l . and the coefficients c K L\mn solve the recursion Eq. 

EZJ, (E23- 

It may be interesting to note that f K Lfi are multiples of derivatives of Legendre 
polynomials (cf. UJ Eqs. 15.2.2, 15.4.4]): 

UM = fT 2 ^ 1 ■ dr l P L+K -i{r). (3.17) 

so that, by Eq. ( 13.161 l, all functions f K L;mn ( r ) are derivatives of the Legendre poly- 
nomials Pi +K _i(r). E.g., for twist 2 (fc = 1), we have 

e lL (p,q,r) = {l + '^(r-l)d r +<L(l+r)d r )P L (r). 
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The next task is to relax the assumption /Xi = jU2 = scalar, and to find and solve 
the analogue of Eq. ( 13. 4t in the general case. This will be necessary in order to com- 
pute the contributions from all insertions of projectors as in Eq. dl.lt by successive 
reduction according to Eq. ( 13.1b . 

For two scalar fields of different dimensions, d\ ^ d%, the ansatz Ex = E k l o 
(xf 2 ) +d2 2 is solved by a scalar polynomial E K i(di,d 2 ,v), homogeneous of de- 
gree 2k + L in di, homogeneous of degree L and harmonic in v, as before, but now 
satisfying the differential equation 

(2(d l Vi)Vi - d l V\ + 2(d 2 V 2 )V 2 - d 2 V 2 2 + {d x - d 2 )(V l - V 2 ))EKL(di,d2,v) = 0. 

Note that the homogeneity conditions require that K is an integer, and that in a GCI 
theory, fields with even twist 2k can arise in the OPE only if d\ — d 2 is even. One 
would therefore have to modify the ansatz when d\ — d 2 is odd. 

Similarly, in the chiral case, the ansatz E d ^' dl = E l ^' dl o {x\ 2 ) dl+dl implies that 
E h u 2 is homogeneous of degree h in <9, and satisfies the differential equation 

(d l V 2 + d 2 V 2 2 + (d l -d 2 )(V l -V 2 ))E h (d u d 2 )=0. 

This is solved by 

# A (ft,4) = £ ^-^+d2) P {P + d x -d 2 ) q 



4 Application: Test of positivity of a 6-point structure 

Recall the positivity problem for the exotic scalar 6-point structures addressed in the 
introduction. We consider here only the simplest example of such a structure, which 
has double poles and is consistent with the constraints due to the requirement that 
the OPE in both the first and last pair of fields starts with twist 2 | fTT| . More general 
double pole structures have been classified in JfQ. 

In im , the leading part of this structure was displayed. In Ifl2l . its "tetraharmonic 
completion" (i.e., the biharmonic completion in both pairs of variables x\,x 2 and 
X5,X(,) was presented in terms of a transcendental function g(s, t). The tetraharmonic 
completion is precisely the twist-2 part in both channels. Unfortunately, however, 
due to a wrong resummation factor, this function g(s,t) was incorrectly computed 
in lfl2l . We shall display the correct function below. 

The leading part of the exotic structure for four scalar fields (j)i,^ 2 , <ps,(j)(, of di- 
mension d and two scalar fields 03 , <j>4 of dimension d' is given by 

( X 15 X 26 X 34 _ 2Xj 5 JC2 3 X4g - 2Xj 5 X2 4 X3g) jj 2 jj 5 6 j 

E(Xl, ■ ■ ■ ,X 6 ) = „ , , , 2 . ( Y 2 \rf'-3 . y 2 Y 2 y 2 r 2 . ( Y 2 W-l ' 

yx n ) x l3 x l4 x 23 x 24 ix 34 j - x 35 A 45- x 36 Jt 46 ^56^ 
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where (•)[<:,/] stands for antisymmetrization. Without loss of generality, we choose 
d = d! = 3. For comparison, we also introduce the following 6-point structure with 
the same symmetries as E, but which has no double poles and appears as part of the 
6-point function of six cubic Wick products of a complex massless scalar free field: 

B{x U --- ,X 6 ) = 2 • ( 2 2 ) 'T'f 7 2 ) ' ( x 2 \2 ' 
\ x 12) \ X U X 23/ [1.2] x 34 \ x -i6 x A5 ' [5,6] ^56' 

The structure B is separately biharmonic in both the 1-2 and 5-6 channels. It turns 
out that the tetraharmonic completion H of B — can be written more compactly 
than that of E given in fl2l . namely 

H(x u ---, X6 ) = (B-| ) ■*(*,»)*(//), (4.2) 



where s = f = and 5' = f' = The condition of biharmonic- 



35 4k_ 

ity amounts to the differential equation [ 12] 

[(1 -tdt){\ +td t +sd s ) - ((l-tdt)+t(2 + td t + sd s ))d s ] g = 

for the function g{s,t). The expansion in a power series in s, g(s,t) = j^g n {t), 
gives the recursion (1 + (n + l)f — t (1 —t)d,)g„ = (1 —td t ){n + td t )g„-\ with go (0 
I . This can be solved in terms of hypergeometric functions, giving 

g(s,t)= ( " + 1 } ' • 2Fi (n,n + 1; 2n + 2; 1 - f). (4.3) 

„ (2n + l)! 

The sum can be performed when the integral representation of the hypergeometric 
functions [TJ Eq. 15.3.1] is inserted, and s,t are expressed in terms of the "chiral 
variables" u± such that s = u + ii- and t = (1 — «+)(l — k_). Then 

g ( Jjt ) = £(n+l)/ / 1 dxx"(l-x)"(l-(l-f)x)-" 

n JO 

'' l \\ — (u + -\- U- — u + iiJ)x 
= I ax 



(1 — u+x){\ — u-x) 

(1 — M+)(l — M-) • log I _" u + +U + — U- — jU 2 + + \u*L 

1 + 2u+U- ■ 



(u+ — M_) 3 



^ L (a + b)((a + b) 2 -l) 



lab 



i.b>l 



(1-^(1 -«_) ^ a - b u> t. (4 .4) 



a,b>0,a+b>0 ° + * 



(In the first line, we corrected a wrong factor of «!, whose presence in [[12] and 
spoiled the subsequent expressions.) 
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Because the twist-2 part is obtained by inserting projections, it must separately 
satisfy Wightman positivity. Of course, we would like to apply the twist-2 inter- 
twiners En °f Sect. [3] in both channels, so that the issue reduces to the positivity 
of tensor-scalar-scalar-tensor 4-point functions. Applying successively the unknown 
intertwiners for the resulting tensor-scalar channels, the problem would be reduced 
to the positivity of the resulting 2-point function, i.e., to the positivity of the numer- 
ical amplitude. 

Since we know the intertwiners E\l, the first step can in principle be done. No- 
tice that it is sufficient to act on the leading part, because it differs from the twist-2 
part by contributions of higher twist, that are annihilated by E\i. Notice also that B 
has the form of a product of two 4-point functions in the variables x\,Xi,x^,x^ and 
in the variables x-^.x^.x^.x^. Therefore, the application of the intertwining differen- 
tial operators in the 1-2 channel and in the 5-6 channel also factorizes. The same, 
however, is not true for E. 

Thus, even the first step at present seems to be too involved to be carried out in 
practice. The second step is at present not possible because we have not yet deter- 
mined the tensor-scalar intertwiners. 

For this reason, we decided to perform only a weaker test of positivity. Namely, 
we restrict the twist-2 structure to 2D, by setting two spatial coordinates to 0. Since 
this essentially amounts to a smaller class of test functions, Wightman positivity 
must still be preserved; but notice that 2D positivity after the restriction is necessary 
but not sufficient to ensure positivity in 4D. 

The intertwining operators in 2D are at our disposal Eq. ( 13.101 . and we have com- 
puted all coefficients (see below). It turns out that the partial wave amplitudes of the 
restricted exotic twist-2 structure B — lE differ from those of the non-exotic struc- 
ture B only by certain signs. This means that E has the same partial wave amplitudes 
as 4-B, except that some of them are absent. 

The non-exotic structure B may itself be indefinite, but we know that it occurs 
in a free-field model, and therefore can be dominated by other positive free-field 
structures, because free fields are manifestly positive. This seems to indicate that the 
restricted exotic structure as well can be dominated by positive free-field structures. 
Thus positivity at the 6-point level alone would not forbid the appearance of this 
structure as part of a 6-point correlation function. 

Let us indicate some details of the actual computations. 

Upon restriction to 2D, u+ and m_ turn into the chiral cross ratios u = X[2 * 34 . 
Moreover, the function B — \E drastically simplifies: 

1 2D 1 (U+ - It-) (u' + - U'_) 

D — —t, ~ 



(x 2 n ) 2 X 2 13 X 2 24 -xl 4 -X 2 35 X 2 46 (x 2 56 f (1 - « + )(l - U_) (1 - U' + )(1 - U'_) ■ 

After multiplication with g(s,t)g(s' ,f'), using Eq. (14.4b . we have 

H 2D 1 v a-b ' a u v a-b b 

( x 2 )2 X 2 2 2 2 2 ( J2 \2 ' a + b + ^ Cl + b + 

\ X 12> x 13^24 x 34 X 3$ X A6\ X 56) a.b>0,a+b>0 a ~ " a,b>0,a+b>0 u T " 
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For the non-exotic structure B, one has instead 



B 



2D 



1 



( x 12) 2x 13 X 24 '*34 ' X 35 X 46( JC 56) 2 a 



u"u b 



b>0,a+b>0 



E 

a,b>0,a+b>0 



u'"u' b . 



Because the sums factorize, the evaluations of the chiral intertwining differential o$- 
emtoisioEi 1± (di c ,di)o(xi ( i±) d in the 1-2 channel^,/ = l,2)and in the 5-6 channel, 
with £/, given by Eq. (13.10b . completely decouple. Actually, because all structures of 



interest are of order xf, 1 d , and therefore only chiral dimensions h > 1 will occur, we 



\d-l 



found it more efficient to work with chiral intertwining operators 1 oD/, ± o (x^ i ±) fl 
where D/,(<9<.,(9/) = (V^ — V/)£fc(dfc,<?/), and adopt a normalization different from 
Eq. EES: 



D h (d 1 ,d 2 ) = 



1 



Thus, we apply i oD h+J ,_ o (x\ 2 ) 2 = I o [D A+ o (x 12 
1 



>£>/>_ °(*i2,-) 2 ]- We find 



■ H 



* 13*24 

W-(l-ft)- 



l 34 



loDh 



l 12 



(xi 3 x 24 



ia+1 



l 34 



(x-x 3 ) h (x-x 4 ) h 



where c a j, = v '"^ — — . Multiplying the two chiral factors and performing the sum 
over a and b gives for the structure B 



loD 



h+,h- 



1 



2 2 i—i 

X 13 X 24 a.b>0,a+b>0 



Cs{h + ,h-) ■ 



4 34,-1 



(x-x 3 ) / ! + (x-x 4 )+ + 



where, by virtue of F(z) :=E B c a jiZ fl = flk-i(l-2z) andft(-l) = (-1) L , 

Cb ft- ) =2Xodd(^), 



(4.5) 



where ^ dd(/ ! ) = 1 if the helicity /; = h + — h- is odd, and zero otherwise. To perform 
the corresponding computation for the sum weighted with as in the structure H, 

one may for b > put G^(z) = fq^fc fl ,/iZ fl , solve the equation zG' + bG = zF' — bF 
by G(z) = F(z) ~2bz~ h fQt b ~ l F(t)dt, and use the orthogonality of the Legendre 
polynomials to conclude G(l) = F(l) = (— 1) /,_1 if h + > One finds 



C H (h+,h-) = sign(ft) • 2^ dd(/j) 



(4.6) 



The same factors arise in the 5-6 channels. Thus, when the 6-point structures B and 
H are reduced in both channels by means of {l x °Dh + ,h- ° t^) 2 ) <S> (ly °Dh' + .ti ° 
(jfjg) 2 ), the result is always a multiple of the same 4-point function 



W/i+A-;/!' ,h'_(x,x 3 ,x 4 ,^) 



4 34,- 



(x - x 3 (x - x 4 )+ + (x 3 - x')+ + (x 4 - 
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The respective coefficients for the structures B and H are 

C B (h+,h-)C B (h' + ,ti_) = 4Xodd(h)Xodd(ti), 

C H (h+,h-)C H (h' + ,ti_) = sign(/i)sign(/i') -^oddMZodd^') , (4.7) 

where /i = h + — /;_, A' = h' + — h'_ are the helicities. 

Because H is the twist-2 part of B — \E, we conclude that (after 2D restriction) all 
partial waves with helicities of equal sign in the 1-2-channel and in the 5-6-channel, 
that are present in B, are absent in the twist-2 part of E, while those with helicities 
of opposite sign arise in the twist-2 part of E with 4 times the coefficient in B. 

It remains to perform the partial wave expansion of the 4-point functions W h+ h _. h i 
Here one may use standard methods, e.g., Ifl5l l6l [TOl . Namely, 

W h+J ,_- h ' + j,'_(x,x 3 ,x 4 ,x) = £ ^+ k h-;h' + ,h'_ ' W ht'hZ;h' + ,h'_^ X3 ' X4 ' X '^ 

k+,k- 

where W k +' k - is the partial wave for the insertion of a projection on the repre 

3 

2 

'h + ,h--Ji',,h'_ 



sentation with scaling dimensions (&+,&_). It turns out that only k± 6 | + No 
contribute. Because of chiral factorization of W h h _. h i y , one has B, 



B h , h , B h _ h , , where the chiral coefficients are determined by the expansion 
1= Y, B k lhl -u n 2 Fi(n + h,n + h';2n + 3;u). 

k= I +n 

The problem of Wightman positivity of the (2D-restricted) structures B and H has 
now been reduced to the positivity of linear combinations of matrices of the form 

P±P &d[B k+ ®B k -]P M P ± , 

where and P± are the projections on the odd resp. positive or negative helicities. 

To be admissible in a QFT, the exotic structure does not need to be separately 
positive, but must only be dominated by other, non-exotic structures that contribute 
to a full 6-point function. Thus, if positivity should fail for H (it certainly does 
for the twist-2 part of E because in this case all diagonal matrix elements vanish), 
one would have to establish a bound for the negative part of the above matrices by 
positive matrices of partial wave amplitudes arising from other structures. 

We have not completed this analysis yet. 

To conclude: the tools are available to test Wightman positivity of 6-point corre- 
lation functions. If a 6-point function involving the exotic structure Eq. ( 14.lt passes 
the test, then it could be a candidate for a nontrivial 4D conformal QFT. 
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